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ABSTRACT 

Durmett  (1955)  developed  a  procedure  comparing  m  treatments  to  one 
control  with  an  exact  overall  type  I  error  of  a  when  all  sampling 
distributions  are  normal.  Sometimes  it  is  desirable  to  compare  m  treatments 
to  controls.  In  particular,  it  is  often  desired  to  compare  m  treatments 
with  two  controls.  For  instance,  several  new  treatments  (e.g.,  pain 
relievers)  could  be  compared  to  two  standard  treatments  (e.g.,  Aspirin  and 
Tylenol).  Alternatively,  a  standard  treatment  could  be  very  expensive, 
difficult  to  apply  and/or  have  bad  side  effects,  making  it  useful  to  compare 
each  new  treatment  to  both  standard  treatment  and  no  treatment  (Placebo) . 

Dunnett's  method  is  expanded  here  to  give  comparisons  of  mean  values  for 

—  *7  -  r 

m  treatments  to  mean  values  for  k^2  controls  at  an  exact  overall  type  I  error 

of  a  when  all  sampling  distributions  are  normal.  Tabled  values  needed  to  make 

v 

exact  simultaneous  comparisons  at  i-.05  are  given  for  k-2.  An  application  is 

made  to  an  illustrative  example  from  the  literature.  C~~  - - 

AMS  SUBJECT  CLASSIFICATION:  62F25,  62F03,  62Q05. 

KEY  WORDS:  Simultaneous  Inference;  Multiple  Treatments;  Multiple  Controls; 
Dunnett's  Procedure. 

I .  INTRODUCTION 


Simultaneous  comparison  of  m  new  treatments  to  one  control  has  been  well 
studied.  The  standard  situation  modeled  is  for  n  observations,  sampled  from 
each  treatment  group  and  the  control  group  with  independence  of  all 
observations,  each  control  group  observation  distributed  N(/rc,  a2),  and  each 

observation  from  treatment  group  j  distributed  N(pt^ ,  a2)  for  j-1 . m. 

Confidence  intervals  for  (p*.  -  #ic)  ,  the  difference  between  each  treatment  mean 

j 

and  t^e  control  mean,  are  desired  (the  treatment  being  declared  different  from 


the  control  if  the  interval  doesn't  cover  zero).  It  is  important  to  limit  the 
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probability  of  one,  or  more,  confidence  intervals  not  containing  the  true 
values  (overall  type  I  error)  to  some  value  a  (usually  or-. 05)  to  prevent 
spurious  findings. 

Often,  the  investigator  is  only  interested  in  a  lower  bound  estimate  for 
the  improvement  each  treatment  gives  over  the  control.  When  larger  means  are 
ueairable,  one  sided  lover  confidence  intervals  of  the  form: 

(jit^-  nc)  >  Lj  (for  Lj  some  function  of  the  data)  j-l,...,m 
can  be  used.  If  smaller  means  are  desirable,  then  one  sided  upper  confidence 
intervals  can  be  used.  From  now  on  WLOG  we  assume  larger  means  are  desirable. 
Sometimes  the  investigator  wishes  to  have  both  upper  and  lower  bounds  for  the 
improvement  each  treatment  gives  over  the  control,  requiring  two  sided 
confidence  intervals  of  the  form: 

Lj  <  (/Jt^  -  /ic)  <  Uj  (for  hj  and  some  functions  of  the  data) 
and  ^jHl,  * . .  ,m. 

Dunnet's  (1955)  well  known  procedure  gives  one  sided  (or  two  sided) 
confidence  intervals  with  a  simultaneous  overall  type  I  error  of  a.  These 
confidence  intervals  have  the  form: 

One  sided:  (jit  -  *ic)  >  (X*  -  X,.)  -  t ™ J2  s/7n  for  j-l,...,m. 

1 

Two  sided:  (/it^  -  jxc)  t  (X*  -  X,.)  ±  t*** Jl  s/Jn  for  j-l,...,m. 

Here  X,.  is  the  sample  average  of  all  control  group  observations,  Xt^ 

is  the  sample  average  of  all  treatment  group  j  observations  j-l,...,m  and  s 
Is  the  pooled  estimate  of  the  standard  deviation. 

Note  that  t(£J  and  t^1  are  constants  that  depend  on  a.  Dunnett's 
confidence  intervals  also  have  the  nice  property  of  symmetry:  the  type  1 


error  for  each  individual  comparison  is  the  same. 

Sometimes  it  may  be  desired  to  simultaneously  compare  m  new  treatments 
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to  k>2  controls.  In  particular,  it  is  often  of  interest  to  compare  a  new 
treatments  to  two  controls .  For  example : 

(A)  -  Two  commonly  used  treatments  (e.g..  Aspirin  and  Tylenol)  are 
each  compared  with  m  new  treatments  (e.g.,  m  new  pain  relievers). 

(B)  -  A  standard  treatment  and  a  placebo  are  each  compared  with  m  new 
treatments.  Two  situations  this  might  arise  in  are:  (i)  The  standard 
treatment  is  expensive,  has  undesirable  side  effects  or  is  otherwise 
impractical  to  use  in  some  situations.  The  new  treatments  do  not  suffer  from 
these  problems  and  thus  under  certain  conditions  will  be  used  if  they  can  be 
shown  to  be  better  than  placebo.  (ii)  The  new  treatments  are  in  a 

developmental  stage.  Promising  treatments  (those  which  work  better  than  no 

/ 

treatment)  may  be  further  developed  and  expanded  testing  will  be  done  of 
treatments  already  performing  better  than  the  standard  treatment. 

Graham,  et  al.  (1988)  compared  a  particular  class  of 
antipyretic/analgesic  treatment  and  a  placebo  to  two  standard  therapies 
(considered  controls  here)  in  a  study  of  undesirable  side  effects.  They 
hypothesized  that  two  standard  treatments  for  cold  symptoms  (Aspirin  and 
Ibuprofen)  had  adverse  effects  on  the  immune  function  (e.g.,  post  virus 
challenge  specific  antibody  response)  and  resulted  in  more  virus  shedding  than 
did  no  treatment  at  all  (placebo).  They  also  hypothesized  that  a  different 
class  of  cold  symptom  treatments,  acetaminophen,  had  less  adverse  effects  on 
the  immune  function  and  caused  less  virus  shedding  than  did  Aspirin  or 
Ibuprofen. 

As  before,  it  is  assumed  that  n  observations  are  sampled  from  each 
treatment  and  control  group  with  independence;  those  from  the  ith  control 
group  (ct)  distributed  N(/ic  ,  a2)  for  i-l,...,k,;  and  those  from  the  jth 
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treatment  group  (tj)  distributed  N(/tt  ,  a2)  for  j— Also  Xt<  is 

J  J 

the  sample  mean  of  all  observations  in  the  jth  treatment  group,  while  Xc  is 
the  sample  mean  of  all  observations  in  the  ith  control  group 

It  is  now  desired  to  give  simultaneous  confidence  intervals  for 
(fir  -  fi.  )  where  j-i",  .  .  .  ,m  and  i— l,...,k.  In  Section  II,  Dunnett's 

Cj  Ci 

methodology  developed  for  k-1  controls  is  expanded  on  to  work  for  k^2 
controls.  The  c  onfidence  intervals  will  be  of  the  form: 


One  sided:  (nt  -  p  )  >  (Xt  -  X  )  -  d(*>y2  s/Jn 
J  c±  cj  ct 

l:  </*t  -  mc,>  <  <xt  -  xCi>  ±  d^j2  s/7H 


(1) 

(2) 


Two-sided: 

where  d***  and  d^are  the  constants  which  make  simultaneous  coverage  of  the 
above  one  and  two-sided  confidence  intervals  equal  to  l-o .  These  confidence 
intervals  also  have  the  property  of  symmetry. 

II.  DEVISING  THE  EXACT  CONFIDENCE  INTERVALS 


The  values  for  d^1  »nd  d1*'  making  the  overall  type  I  error  of  the  above 
one  and  two-sided  confidence  intervals  equal  to  a  are  derived  using  order 
statistics  of  the  standardized  control  means  along  with  mutual  independence  of 
control  means,  treatment  means  and  s.  Two  sided  confidence  interval  results 
are  more  difficult  to  obtain  and  are  dealt  with  first, 
a.  Two  Sided  Confidence  Intervals 

Define  the  standardized  treatment  group  means  and  control  groups  means 


as : 

T,  -  Oi  -  xt  )  fr/o  for  j  -  l...,m  (3) 

J  J 

ci  -  (Hc  ~  Xc  )  yH/a  for  i  -  l...,k.  (4) 

ci  ci 

Now  all  two-sided  confidence  intervals  given  in  (2)  contain  the  true 
differences  if  and  only  if: 
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|Tj  -  CA  |  <  d<2,72  •  s/a  for  all  i  and  j.  (5) 

If  W  and  V  are  the  maximum  and  minimum  values  among  C1,C2,  .  .  .  ,Ck  then  (5)  is 
true  if  and  only  if : 

T,  -  W  >  -  •  s/a  for  all  j-1 . m  (6) 

and 

Tj  -  V  <  d(*V2  •  s/a  for  all  j-l,...,m.  (7) 

Now  Tlt  T2,  .  .  .  ,T.  and  Clf  C2,...,Ck  are  all  independent  and  identically 
distributed  N(0,1)  variables.  The  joint  density  of  the  order  statistics  (W, 

V)  is  easily  obtained  from  a  well  known  formula,  e.g.  in  Mood,  Graybill  and 
Boes  (1974)  and  is: 

f(V— v,  W-w)  -  2  ^(v)  ^(w)  [$(w)  -  4>(v)  ]k~2  I  { w> v ) 
where  9>  and  <$>  are  the  standard  normal  cumulative  distribution  and  standard 
normal  density  function  respectively.  If  S0-s/cr  and  r  is  the  degrees  of 
freedom  of  s,  then  r(Sc)2  has  a  X2r  distribution.  Also  (W,  V)  is  independent 
of  Sc  and  of  Ti,...,Ta.  For  any  fixed  values  S0-s0,  V-v  and  W-w,  (w>v)  the 
probability  that  {  |  Tx  -  CA  j  <  J2  S„  d(2>  )  for  all  i  -  1,  .  .  .  ,k  is  the  same 
as  the  probability  that  {  -  v  <  Jl  sc  d(2)and  Tx  -  w  >  J2  sc  d(2)  )  .  This 

probability  equals: 

Max  [0,  ($(v  +  Jl  scd<2>)  -  *  (w  -  72  sc  d(2>  )  ]  .  (8) 

By  independence  of  the  TJ:j-l,...,m  and  still  conditioning  on  W-w,  V-v 
and  S0-s0,  the  probability  that  {  Ts  -  v  <  72  sc  d(2)  and  Tx  -  w  >  -  Jl  se  d(2)  } 
for  all  j  is: 

[Max  [0,  <*  (v  +  72  sQ  d<2>)  -  *(w  -  Jl  sc  d<2> ))]]-. 

Using  the  independence  of  (W,  V)  from  Se  and  integrating  over  their 
joint  density  gives  the  unconditional  probability  of  coverage  by  the  intervals 
in  (2)  as: 
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/-I  J-«  A<rso2)  *(v)  *(w)  [4>(w)  -  4>(v)]k"2 
[Max  [0,  4>(v  +  Jl  d(2>  sc  )  -  4>(w  -  Jl  dt2>  s0)]]"  d(rs02)  dv  dw.  (9) 
where  \frr(*)  is  the  standard  X2r  density  function.  This  three  dimensional 
interval  can  be  evaluated  using  a  powerful  computer  such  as  a  Cray,  FPS  264  or 
an  IBM  3030. 

When  the  probability  in  (9)  equals  (1-a) ,  then  the  probability  that  all 
the  intervals  in  (2)  simultaneously  contain  the  true  differences  is  (1-a) . 

This  means  that  the  overall  Type  I  error  is  a.  By  iteration,  from  initial 
guesses,  the  value  of  d(2)  which  will  set  the  probability  in  (9)  equal  to 
(1-a)  can  be  found. 

Note  that  when  k  is  two  (i.e.,  there  are  two  control  groups)  then 
($(w)  -  $(v)]k-2  -  1,  and  so  this  term  drops  out  of  formula  (9).  When  a  is 

known,  or  asymptotically  when  r  is  large,  Se  can  be  set  to  1  eliminating  the 
inner  integrand  and  reducing  (9)  to  a  two  dimensional  integral  that  can  be 
easily  evaluated  on  most  mainframe  computers.  There  do  not  appear  to  be  any 
asymptotic  simplifications  of  (9)  as  m  and/or  k  get  large. 

b.  One  Sided  Confidence  Limits 

The  one  sided  lower  confidence  limits  given  in  (1)  will  all  contain  the 
true  values  if  and  only  if 

(Tj-Ci)  >  Jl  S0  d<2)  for  all  j-1 . m  and  i-1 . k.  (10) 

which  is  true  if  and  only  if 

(Tj— V)  >  Jl  SD  d<”  for  all  j  -  l,...,m.  (11) 

With  the  same  procedure  used  for  2  sided  confidence  intervals,  the 
simultaneous  coverage  probability  of  the  one  sided  upper  limits  given  in 
(10)  and  (11)  can  be  shown  to  equal  the  following  integral: 

LI  1 1  tft(rs02)  rf(v)  [1  -  *(v)]k“l  [*(v+/2d<i>  so)]»  d(rs02)  dv.  (12) 
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This  two  dimensional  integral  is  easily  evaluated  on  most  mainframe  computers. 
The  value  of  d^ which  makes  the  overall  Type  I  error  equal  to  (1-a)  is 
obtainable  through  iteration  from  initial  guesses. 

When  a  is  known,  01  asymptotically  when  r  is  large,  then  SG  can  be  set 
to  1  eliminating  the  inner  integrand  and  reducing  (12)  to  a  one  dimensional 
integral.  Note,  that  by  symmetry,  values  for  doused  for  one  sided  upper 
confidence  limits  will  be  the  same  as  those  needed  for  one  sided  lower 
confidence  limits. 

HI.  TABLED  VALUES 

Tables  1  and  2  give  values  of  d(£5  and  ,  respectively  for  k— 2 ,  a— .05, 

m-1,2 . 10,  and  various  values  of  r.  Values  for  k  >  3  could  be  just  as 

easily ^ calculated  but  are  not  given  here  since  it  is  felt  that  use  of  3  or 
more  controls  will  not  be  common. 

The  tables  were  constructed  by  numerical  evaluation  of  the  intergrals  in 
(9)  and  (12)  to  an  accuracy  of  0.0001.  The  integration  was  performed  with  a 
Fortran  program  using  an  adaptive  Romberg  algorithm  that  was  taken  from 
Rabinowitz  (198a),  aiiu  altered  to  greatly  reduce  the  number  of  computational 
steps.  Secant  iteration  was  used  to  obtain  d(£} and  d‘5*  to  the  nearest  fourth 
decimal  place.  The  results  were  then  rounded  to  the  highest  third  decimal 
place  and  included  in  the  table.  Good  approximations  of  values  of 
dm 

(or  d^5)  associated  with  r  not  in  the  tables  are  obtainable  through 
standard  linear  interpolation  of  d^*  (or  d1^)  associated  with  tabled  values 


of  r. 
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The  improvement  from  using  the  exact  value  for  k-2  over  using  the  best 
available  upper  bound  (Tukey  or  Bonferroni)  depends  on  a,  r  and  m.  Generally, 
the  improvement  increases  with  m  and  decreases  with  r.  For  instance,  the 
exact  value  for  d(£5  with  (m-1,  r—*>)  is  1.917  which  is  2.19%  smaller  than  is 
the  Bonferroni  upper  bound  of  1.960.  When  (m-1,  r-5)  the  exact  value  for  d(^} 
is  2.441,  which  is  5.06%  smaller  than  is  the  Bonferroni  upper  bound  of  2.571. 
When  (m-5,  r— »)  the  exact  value  is  2.487,  which  is  3.45%  smaller  than  is  the 
Bonferroni  upper  bound  of  2.576.  Finally,  when  (m-5, r-5)  the  exact  value  of 
3.469  is  13.96%  smaller  than  is  the  Bonferroni  upper  bound  of  4.032. 

For  small  a,  the  improvement  decreases  as  a  decreases,  due  to  the 
phenomenon  of  the  Bonferroni  approximation  converging  to  the  exact  value  as  a 
decreases  which  was  noted  by  Dunn  (1958).  Also,  as  k  increases,  the 
improvement  of  the  exact  value  over  upper  bounds  will  increase  for  the  same 
reason  it  does  as  m  increases. 

IV.  ILLUSTRATIVE  EXAMPLE 

Below  is  an  illustrative  example  of  two  treatments  being  simultaneously 
compared  to  a  control  and  a  placebo.  It  is  adapted  from  one  given  by  Villars 
(If 51)  and  used  in  Dunnett  (1955).  The  data  represent  measurements  on  tne 
breaking  strength  of  untreated  fabric  (placebo) ,  fabric  treated  by  an 
expensive  standard  method  and  fabric  treated  by  two  proposed  less  expensive 


methods . 
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Breaking  Strength  (lbs.) 

Ci-No 

Trt 

C2— Std  Trt 

Tj— New  Trt  1 

T2-New  Trt  2 

50 

55 

55 

55 

Observations 

41 

64 

49 

47 

61 

£2 

48 

Mean 

45 

61 

52 

50 

Variance 

01 
«-  A 

27 

9 

14 

Here  ro-2  and  n-3.  The  Pooled  variance  estimate  is  s2  -  19  and  the  degrees  of 
freedom  is  8.  The  estimated  standard  error  of  a  difference  between  two  means 
is  s  fT/n  -  3.56.  Symmetric  two  sided  confidence  intervals  for 
(/i,.  -  p  )  with  an  overall  Type  I  error  of  a  would  take  the  form  given  by 

j  i 

(2): 


(X,.  -  X  )  ±  d<2,72/n  s  for  all  i  and  j. 

j  Ci 

Symmetric  one  sided  lower  confidence  limits  would  be  given  by  (1): 
(Xr  -  X,,  )  -  d(*J  JT/x\  s  for  all  i  and  j  . 

j  i 


The  exact  values  for  dl2)  and  dll)  obtained  from  Tables  2  and  1  along 

.05  .05 

with  commonly  used  upper  bounds  are  given  below. 

Upper  Bounds 


Exact 

( Bonferroni 

Tukey 

Schef fe ] 

d<2> 

3.053 

3.210 

3.200 

3.492 

.05 

d<n 

.05 

2.588 

2.751 

N.  A. 

N.A. 

Ninety-five  per  cent  simultaneous  confidence  intervals  for  (treatment-control) 
differences  derived  from  the  exact  and  best  upper  bounds  are  now  given. 
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Confidence  Intervals  C 95% 'i 

Comparison 

Two 

Sided 

One  Sided 

(Lower) 

Exact 

Tukey 

Exact 

Bonferroni 

'"'l 

(-3.87,  17.87) 

(-4.39,  18.39) 

>  -2.22 

>  -2.96 

Mr 

- 

(-19.87,  1.87) 

(-20 . 39 ;  2.39) 

>  -18.22 

>  -18.96 

“1 

c2 

M*2 

"^1 

(-5.87,  15.87) 

(-6.39,  16.39) 

>  -4.22 

>  -4.96 

"*2 

(-21.87,-0  13) 

(-22.39,  0.39) 

>  -20.22 

>  -20.96 

The  exact  method  finds  che  second  new  treatment  to  be  significantly 
worse  than  the  standard  treatment  when  making  two  sided  comparisons  at  a-. 05. 
This  statistical  difference  is  not  seen  at  o-.05  using  the  closest  upper 
bound. 

V.  VARIATIONS  ON  ASSUMPTIONS 

Formulas  (9)  and  (12)  can  be  altered  to  allow  for  unequal  sample  sizes 
and  different  variances  within  the  control  and  treatment  groups.  One  such 
variation  would  be  having  observations  sampled  for  each  treatment  group  and 
observations  sampled  for  each  control  group  with  control  and  treatment 
variances  the  same.  Under  these  assumptions,  formula  (9)  becomes: 

L«  /’  V’t (rs02)  2  <Hv)  4>( w)  [4>(w)  -  $(v)]k~2  [Max  [0,  (4>(v  +  ( (nc+nt)/nt)h 
d(2)s0)  -  $(w  -  ((r^+nj/njh  d(2)s„))]“  d(rsc2)  dv  dw. 

When  given  a  total  of  n„  observations  to  be  divided  among  m  treatment 
groups  and  k  control  groups.  It  can  be  shown  by  calculus  (as  Dunnett  (1955) 
did  for  k-1)  that  the  optimal  allocation  of  a  fixed  number  of  observations 
(n„)  to  minimize  the  variance  of  (X,.  -  X  )  is:  ti*  -  7k  i\,  /  ((7k  +  7m)  m) 

j  i 

and  nc  -  7“  ^  /((7k  +  To)  k) . 
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VI.  VARIATION  ON  THE  CONFIDENCE  INTERVALS 


Comparing  Both  Control  Groups  to  Each  Other 

In  case  (A)  of  the  introduction,  it  may  sometimes  be  of  interest  to  also 

symmetrically  compare  the  two  control  groups  to  each  other,  as  well  as, 

comparing  each  treatment  to  each  control,  with  an  overall  Type  I  error  of  a. 

To  do  this  requires  building  a  two  sided  confidence  interval  for  (u_  -  u_  ) 

C1  c2 

of  the  form  (X  —  X  )  ±  d_  Jls  in  addition  to  the  one  or  two  sided 
C1  c2 

confidence  intervals  for  the  difference  in  (1).  If  the  intervals  for  the 
differences  in  (1)  are  two  sided,  then  the  formula  for  simultaneous  coverage 
is : 


/-I  /;  /  g  ^r(rso2)  2  ^(v)  4(v)  (*<v  +  -Z2  d(2)s0)  -  $(*'  -ft  d(2)sj)* 
dv  d(rs02)  dw 

where  e  -  w  -  d<2)./2  s0  If  the  intervals  for  the  differences  in  (1)  are  one 
sided,  then  the  formula  for  simultaneous  coverage  is: 

Jo  J”  lMr s02)  2  ^(v)  ^(w)  ($(v  +  72  d^’s,,)*  dv  d(rsc2)  dw 
where  e  -  w  -  d^  ft  s0 
Other  Variations 

It  may  be  desired  to  have  two  sided  confidence  intervals  for  the 
differences  between  all  treatments  and  the  first  control  and  one  sided 
confidence  limits  for  the  differences  between  all  treatments  and  the  second 
control.  Sometimes  the  mean  effect  of  the  first  control  may  be  known  while 
that  of  the  second  control  is  not.  Both  of  these  previous  situations  may 
occur  when  the  second  control  is  a  placebo.  For  each  of  these  situations  it 
is  possible  to  produce  symmetric  confidence  intervals  having  an  overall  Type  I 
error  of  a  using  modif ications  of  the  previous  methodology.  The  formulas 
needed,  however,  are  quite  complicated. 
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TABLE  1:  Exact  Values  of  for  k  -  2  Controls,  m  New  Treatment  and  r  Degrees  of  Freedom. 


/ 


vO 

00 

o 

r-. ' 

CM 

vi) 

H 

I/O 

m 

ON 

o 

m 

H 

OO 

VO 

in 

m 

m 

VO 

o 

r** 

CO 

00 

\o 

00 

«H 1 

VD 

CM 

00 

uo 

CM 

o 

00 

n- 

m 

<* 

On 

in 

H 

CO 

ON 

rH 

00 

Vi) 

CO 

CM 

CM 

fH 

rH 

o 

o 

O 

o 

o> 

ON 

ON 

ON 

OO 

00 

OO 

r*v 

r-* 

VO 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

VO 

ON 

ON 

CM 

H 

o 

in 

in 

o 

00 

ON 

CM 

OO 

<r 

CO 

CM 

iH 

CM 

CO 

m 

H 

OO 

CM 

rH 

CO 

CM 

OO 

H 

ON 

VO 

CO 

♦H 

o 

VO 

CM 

00 

3 

o 

VO 

ON 

oo 

m 

•O 

co 

CM 

H 

rH 

o 

© 

o 

On 

On 

ON 

ON 

On 

ON 

OO 

CO 

r- 

rv 

vO 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

co 

co 

co 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

oo 

CO 

CM 

in 

in 

CO 

CO 

CO 

ON 

o 

vO 

in 

vO 

o 

VO 

CO 

CM 

CM 

CO 

n 

oo 

H 

3 

CM 

vO 

oo 

CM 

C*”N 

ON 

3 

CM 

o 

3 

VO 

CM 

00 

o 

vO 

CO 

00 

in 

CO 

CM 

rH 

H 

o 

o 

ON 

On 

ON 

ON 

ON 

OO 

oo 

00 

00 

r- 

vO 

vO 

co 

co 

CO 

co 

CO 

CO 

CO 

co 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

o 

CM 

3 

VO 

vO 

oo 

ON 

VO 

On 

n 

in 

f—f 

in 

in 

r^. 

ON 

CO 

VO 

f — 1 

in 

o 

o 

CM 

vO 

rH 

r^* 

3 

rH 

ON 

r-* 

in 

CM 

«— i 

co 

On 

VO 

CM 

ON 

vO 

3 

CO 

CM 

r—l 

o 

o 

ON 

ON 

On 

oo 

oo 

00 

00 

O0 

oo 

n- 

vO 

VO 

VO 

in 

co 

CO 

CO 

CO 

co 

CO 

CO 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

ON 

o 

I— 1 

CM 

o 

vo 

in 

oo 

VO 

VO 

ON 

-3- 

o 

00 

rH 

00 

CO 

00 

*3- 

CO 

CO 

in 

o 

in 

i— 1 

00 

in 

CO 

rH 

On 

oo 

m 

CM 

00 

iH 

vO 

n 

CO 

CM 

»— l 

o 

o 

ON 

On 

00 

oo 

oo 

00 

r-. 

n- 

r- 

vo 

vO 

vO 

n 

in 

-* 

co 

CO 

CO 

CO 

co 

CO 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

On 

rH 

ON 

m 

>3 

o 

VO 

o 

H 

VO 

•3 

in 

ON 

in 

CM 

rH 

m 

o 

vO 

CM 

ON 

SO 

CO 

•3 

f*-* 

CM 

r- 

-3 

T-l 

CO 

vO 

3 

CM 

fH 

o 

ON 

m 

CM 

in 

in 

i-M 

00 

in 

-3 

CM 

f— 1 

O 

ON 

ON 

00 

00 

O0 

r>. 

VO 

vo 

vO 

oo 

in 

m 

<3 

CO 

CO 

CO 

co 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CO 

00 

in 

n* 

f-H 

o 

ON 

vo 

00 

in 

n 

CM 

n* 

CM 

vO 

CO 

o 

00 

VO 

m 

CO 

3 

CM 

CO 

CM 

n- 

3 

rH 

ON 

in 

3 

m 

3 

o 

r- 

3 

o 

3 

f-H 

O’ 

CO 

f-H 

o 

ON 

oo 

OO 

r*^ 

n- 

vo 

vo 

vo 

vO 

vo 

VO 

vo 

m 

in 

in 

3 

3 

3 

CO 

co 

co 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

3 

VO 

3 

m 

in 

00 

rH 

f-H 

VO 

3 

VO 

o 

vo 

3 

co 

3 

CO 

CO 

3 

VO 

00 

o 

in 

00 

r- 

On 

CO 

00 

in 

CM 

ON 

in 

3 

CM 

i-H 

o 

ON 

VO 

CO 

O 

3 

CM 

CO 

H 

ON 

oo 

r** 

VO 

vo 

VO 

in 

in 

in 

m 

m 

m 

in 

3 

3 

3 

3 

co 

CO 

CO 

cO 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

VO 

CM 

VO 

00 

VO 

vO 

3 

00 

r*- 

CO 

CM 

ON 

vO 

ON 

CO 

CM 

CO 

ON 

ON 

in 

in 

O0 

CO 

ON 

VO 

CO 

f-H 

ON 

00 

vO 

m 

3 

co 

co 

o 

m 

CM 

o 

CM 

00 

vo 

in 

in 

3 

3 

3 

3 

CO 

CO 

CO 

CO 

CO 

co 

CM 

CO 

CM 

CM 

CM 

CM 

f-H 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

rH 

r*« 

00 

00 

o 

rH 

OO 

ON 

CO 

ON 

r*N 

n- 

3 

H 

3 

00 

OO 

© 

rH 

CO 

m 

3 

co 

vO 

fH 

00 

in 

CM 

o 

On 

vO 

in 

3 

3 

CO 

CM 

o 

ON 

n* 

in 

cO 

?H 

3 

CO 

CM 

CM 

rH 

rH 

rH 

rH 

o 

o 

o 

o 

O 

o 

o 

o 

o 

ON 

ON 

ON 

On 

On 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

rH 

rH 

iH 

fH 

»H 

OfHcMco<rinvor-.oooNO<*ooc>cM 

U  I0v0r^  .'*0>t— irHrHrHfHiHrHrHrHfHCMCMCO^VOfHg 


TABLE  2:  Exact  Values  of  d^^  for  k  -  2  Controls,  ra  New  Treatments  and  4  Degrees  of  Freedom. 
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